Comparison between numerical relativity and a new class of post-Newtonian 
gravitational-wave phase evolutions: the non-spinning equal-mass case 
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We compare the phase evolution of equal-mass nonspinning black-hole binaries from numerical rel- 
ativity (NR) simulations with post-Newtonian (PN) results obtained from three PN approximants: 
the TaylorTl and T4 approximants, for which NR-PN comparisons have already been performed in 
the literature, and the recently proposed approximant TaylorEt. The accumulated phase disagree- 
ment between NR and PN results over the frequency range Mui — 0.0455 to Mui = 0.1 is greater 
for TaylorEt than either Tl or T4, but has the attractive property of decreasing monotonically as 
the PN order is increased. 
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I. INTRODUCTION 

The current interferometric gravitational-wave detec- 
tors [l|, Q, El have reached design sensitivity, and have 
finished taking data in the S5 science run. Gravitational 
waves from coalescing black-hole binaries will be among 
the strongest that one hopes to find in the detector data, 
and data analysts are searching for them by performing 
matched filtering against template banks of theoretical 
waveforms. A broad class of theoretical templates can 
be produced using post-Newtonian (PN) approximation 
techniques, which are expected to be valid during a bi- 
nary's slow inspiral. 

There exist several prescriptions to compute 
gravitational-wave (GW) templates at different or- 
ders of a PN expansion. Each prescription, termed a 
PN approximant, provides a slightly different GW phase 
and hence frequency evolution. For the purpose of GW 
data analysis, it is important to know which approach, 
at which PN order, best approximates the true phase 
evolution of GW signals from astrophysical black-hole 
binaries. If PN expansions had simple convergence 
properties and could be taken to arbitrarily high order, 
one would expect all approaches to converge to the same 
(presumably correct) result. However, PN calculations 
do not currently go beyond 3.5PN order, and the 
convergence properties are far from clear. By 3.5PN 
order, we mean all the corrections up to the relative 
order x 7 / 2 , where x = (Mujb) 2 ^ 3 , M and cj& are the 
binary's total mass and orbital angular frequency. If we 
consider the number of GW cycles in a given frequency 
window, for example, the number usually docs not 
change monotonically as the PN order is increased. 

One way to check the physical accuracy of the PN re- 
sults is to compare with fully gencral-rclativistic numer- 
ical simulations over the last several orbits of a binary's 
evolution. Recent breakthroughs in Numerical Relativ- 
ity (NR), reported in Ref. [I Ell], 

have made it possible 
to simulate many orbits before merger @, H, , with the 
lar gest number of orbits currently achieved being about 
15 . These simulations allow a detailed comparison 



with various PN approximate GW phase evolutions. Re- 
cent studies have shown that standard PN approximants 
give good phase agreement with NR results up to a few 
orbits before merger 0, 0, El; for example, for a non- 
spinning equal-mass binary, the accumulated phase dis- 
agreement between NR and PN results that use the Tay- 
lorTl approximant, for the seven orbits before the GW 
frequency reaches Mui = 0.1 (about two orbits before 
merger), is less than 1 radian. 

A number of PN GW template families have been pro- 
posed for use in GW searches (see, for example, jll|). 
and have been implemented in the LSC Algorithms Li- 
brary (LAL) [13] for that purpose. We focus here on 
those based on the Taylor approximants. These approxi- 
mants model GWs from compact binaries inspiraling due 
to radiation reaction, and the phase evolution equations 
are expanded in terms of x. Recently one of us proposed 
a new class of templates [lU in which the phase evo- 
lution d(f>/dt is now given by a PN expansion in terms 
of the binding energy (see Section HTT1 for more details). 
Following the LAL terminology, we refer to this new ap- 
proximant as TaylorEt. The TaylorEt template has the 
attractive feature that the number of accumulated GW 
cycles changes monotonically as we increase the reactive 
PN order, which is not true for the other Taylor approxi- 
mants. Further, TaylorEt templates have the same com- 
putational cost as, for example, the TaylorTl and T4 
approximants. Additionally, it has been recently demon- 
strated, while restricting radiation reaction to dominant 
quadrupole contributions, that the TaylorEt templates 
are more efficient than TaylorTl and T4 in capturing 
GWs from inspiraling compact binaries having orbital 
eccentricity [lj]. Therefore, it is interesting to compare 
the GW phase evolution under the TaylorEt prescription 
at various reactive PN orders with the GW phase evolu- 
tion in numerical simulations. This is what we pursue in 
this paper. 

In this investigation, we compare the GW phase evolu- 
tion in numerical simulations of equal-mass non-spinning 
black-hole binaries that last about nine orbits before 
merger with its counterparts obtainable from TaylorTl, 
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T4 and Et prescriptions at various PN orders. It is im- 
portant to emphasize that our NR simulations have low 
initial orbital eccentricity e ~ 0.0016 and an accumu- 
lated numerical uncertainty of less than 0.25 radians in 
GW phase evolution. We observe that the accumulated 
GW phase difference between TaylorTl, T4 and NR os- 
cillates as we increase the reactive PN order from 2PN to 
3.5PN in steps of 0.5PN order. However, differences in 
accumulated GW phases between TaylorEt and NR de- 
scriptions decrease monotonically as we increase the reac- 
tive phase evolution from 2PN to 3.5PN order in steps of 
0.5PN order. This implies that the least GW phase dif- 
ference associated with TaylorEt occurs at 3.5PN order 
and it is A<f> = 1.18 radians. The observed convergence 
of the TaylorEt GW phase evolutions towards their NR 
counterpart and the tolerable accumulated GW phase 
difference for TaylorEt at 3.5PN order are the two main 
conclusions of the present paper. 

In Sectionjlllwe briefly describe the numerical methods 
used to produce the NR waveforms. Section IHII details 
the construction of different PN prescriptions for the GW 
phase evolution of the inspiral of an equal-mass binary 
modeled by nonspinning point-particles. Section [TV] ex- 
plains how we make contact between NR and various 
PN descriptions for GW phase evolutions, and we dis- 
play our results and provide explanations for our obser- 
vations. Concluding remarks and future directions are 
given in Section fVl 



II. NUMERICAL METHODS AND 
WAVEFORMS 

The numerical simulations were performed with the 
BAM code HE EH, where fourth-order accurate deriva- 
tive operators were replaced by sixth-order accurate spa- 
tial derivative operators in the bulk, as described in [8|. 
The code started with black-hole binary puncture ini- 
tial data [13, EH generated using a pseudo-spectral code 
[T^ | , and evolved them with the variant of the moving- 
puncture dEHl version of the BSSN [HH3] formulation 
of the 3+1 Einstein evolution equations. The gravita- 
tional waves emitted by the binary were calculated from 
the Newman-Penrose scalar \l/4, using the implementa- 
tion described in fTsT ]. 

The set of simulations that we discuss here consists 
of binaries that begin with a coordinate separation of 
D = 12M. The initial momenta are chosen according 
the prescription given in [HJ , which lead to inspiral with 
a minimal eccentricity of e < 0.0016. For comparison, 
a second set of simulations uses initial momenta that 
lead to a larger eccentricity of e ss 0.008. Simulations 
were performed at three resolutions. The results were 
seen to converge consistent with sixth-order accuracy and 
were subsequently Richardson extrapolated with respect 
to resolution. Gravitational waves were extracted at five 
extraction radii, and the GW amplitude was extrapolated 
with respect to extraction radius to estimate the GW am- 



plitude as measured at infinity. We use the GW phase as 
measured at the largest extraction radius (R ex = 90 M). 
This procedure is described in detail in 0, where the 
results of these simulations were first reported. 

The resulting gravitational waveform has an uncer- 
tainty in the amplitude of less than 2%, and an accumu- 
lated uncertainty in the phase of less than 0.25 radians. 
In the following sections we will focus on comparing the 
GW phase with that calculated by PN methods, which 
we will now summarize. 



III. VARIOUS PRESCRIPTIONS FOR GW 
PHASE EVOLUTIONS IN PN RELATIVITY 

During the inspiral of a compact binary, the associated 
temporal GW phase evolution can be accurately mod- 
eled using the PN approximation to General Relativity. 
In this approximation, it is customary to consider inspi- 
raling astrophysical compact binaries undergoing adia- 
batic inspiral along circular orbits due to the emission of 
gravitational waves. In recent years a number of com- 
putational efforts provided four particularly valuable PN 
expressions that are essential for GW astronomy with 
inspiraling non-spinning astrophysical compact binaries. 
These four quantities are the 3PN accurate dynamical 
(orbital) energy £(x), expressed as a PN series in terms 
x, the 3.5PN accurate expression for GW energy luminos- 
ity C(x) and the 2.5PN amplitude corrected expressions 
for h + (t) and h x (t) , written in terms of the orbital phase 
<f> and x [IE SE [23, SI]. 

GW data analysis groups focusing on inspiraling com- 
pact binaries employ these inputs to construct various 
types of search template. In this paper, as mentioned 
earlier, we will first consider the TaylorTl and T4 ap- 
proximants. These two template families employ the fol- 
lowing expression for the restricted PN waveform 



h(t) oc x(t) 2/3 cos2 0(t) . 



(1) 



where the proportionality constant may be set to unity 
for non-spinning compact binaries. At a given PN order, 
the two families provide two slightly different ways to 
compute x(t) and <j>(t). To obtain the GW phase evolu- 
tion <j)(t) in the TaylorTl approximant, one numerically 
solves the following two differential equations: 



dtj>(t) 

dt 
dx(t) 

dt 



u b (t) = 
-C(x) 



GM 



r 3/2 



(2a) 
(2b) 



and in this section, we do not employ geometrized units. 
This implies that to construct TaylorTl 3.5PN search 
templates, one needs to use 3.5PN accurate C(x) and 
3PN accurate £(x), respectively The explicit expressions 
for these quantities, in the case of equal mass compact 



3 



binaries, read 



solving the following coupled differential equations. 
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where 7 is the Euler gamma. 

An alternative PN approximant, TaylorT4, has re- 
cently been introduced [lOj. The TaylorT4 approximant 
is obtained by Taylor expanding the right hand side of 
Eq. (|2b[) for dx/dt and truncating it at the appropri- 
ate reactive PN order. Therefore, to construct TaylorT4 
3.5PN search templates, the following set of differential 
equations are numerically integrated: 



dcj){t) 

dt 
dx{i) 

dt 
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The Taylor Tl approximant has been c omp ared with sev- 
eral sets of numerical simulations lldl | and found to 
agree with the NR phase to within one radian over the 
GW frequency range Mcu = 0.05 and Mcu = 0.1, when 
matched at Mcu = 0.1. (Note that the GW frequency 
is related to the binary's orbital frequency by roughly 
cu = 2u>b.) The Taylor T4 approximant has been found 
to a gree within 0.06 radians when compared in the same 
way [1 01 ] . 

Let us now describe the TaylorEt approximant. The 
restricted PN waveform reads 



h(t) cx £{t) cos2 0(i) . 



(5) 



The time evolution for <j)(t) and £ (t) are obtained by 
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where £ = — 2 c 2 and [i being the usual reduced mass. 
In this paper we keep defi/dt to 3PN accuracy (and this 
is the highest PN order available for Eq. (|6"a|) associated 
with compact binaries in PN accurate circular orbits). 
The values of £ corresponding to any initial and final GW 
frequencies can be numerically evaluated using the right- 
hand side of Eq. (|6"a|) . This is possible due to the fact that 
for GWs from compact binaries, having negligible orbital 
eccentricities, the frequency of the dominant harmonic is 
/gw = w&/2. 

At first sight the only difference between TaylorTl/T4 
and Et is a different choice of expansion variable. This 
in itself is interesting to consider: we would like to ver- 
ify that different (valid) choices of expansion variable do 
not dramatically change the predictions from PN approx- 
imants. In addition, however, the TaylorEt approximant 
can be easily generalized to eccentric binaries, and in- 
deed the lower order terms in Eq. (|6ap are responsible 
for precession of an eccentric binary. 

It should be noted that GW phase and frequency evo- 
lutions are prescribed in certain parametric PN accurate 
ways in the TaylorEt approximant. This is achieved in 
the TaylorEt approach by prescribing d<j)/dt, governing 
the PN accurate conservative orbital phase evolution, in 
terms of £ and then numerically imposing secular changes 
in £ , due to the emission of GWs, with the help of far- 
zone GW energy flux, expressed in terms of £. 

In the next section, we compare GW phase evolutions 
predicted by NR simulations and TaylorTl, T4 and Et 
approximants at various PN orders in a given GW fre- 
quency window. 



IV. MAKING CONTACT BETWEEN NR AND 
PN GW PHASE EVOLUTIONS 

In FigureHJ we plot GW phase differences, with respect 
to NR simulations in a given GW frequency window, at 
four post-Newtonian orders, namely 2PN, 2.5PN, 3PN 
and 3.5PN, associated with the three different PN ap- 
proximants described earlier. The GW frequency window 
we employ is between Mcu = 0.0455 and Mcu = 0.1, and 



4 



we line up the NR and PN GW phases and frequencies 
at Muj = 0.1. The top panel refers to TaylorTl and we 
observe the usual fluctuating differences in GW phases 
with respect to the PN order. The center panel shows a 
comparison with the Taylor T4 approximant. The good 
agreement at 3.5PN order that was reported in [H| is 
clearly visible. Again, we observe that the PN-NR phase 
disagreement fluctuates with respect to PN order. Fur- 
ther, substantially different GW phase evolutions in com- 
parison with their NR counterpart for TaylorTl and T4 
at 2.5PN order are also observed. 

In the bottom panel, we focus our attention on the 
TaylorEt approximant. When we increase PN accuracy 
of the reactive dynamics, given by Eq. (|6b[) . from 2PN 
to 3.5PN in steps of 0.5PN, the GW phase difference, 
i.e. 0pn — <^NR; decreases in a monotonic manner. This 
implies that 0pn — </>nr has the lowest value at the 3.5PN 
order, and is given by A(j> = —1.18 radians. This GW 
phase difference at 3.5PN order is definitely more than 
its counterparts arising from TaylorTl and T4 (Acp = 
0.6 and 0.05 respectively). However, it is intriguing that 
the TaylorEt approximant seems to display monotonic 
convergence towards the NR phase. We have checked 
that monotonic convergence to NR is independent of the 
choice of our matching frequency. 

It is possible to explain, using PN based arguments, 
the comparatively large gap one finds between 2.5 and 
3PN orders in the bottom panel of Fig. [T] Recall that 
the conservative phase evolution in the TaylorEt approxi- 
mant is 3PN accurate (in other words, we always use 3PN 
accurate expression for dcj)/dt, given by Eq. (|6"a|) ) and 
this requires a complete knowledge about the 3PN ac- 
curate conservative orbital dynamics. Further, it should 
be clear that the reactive GW phase evolution in Tay- 
lorEt is governed by Eq. (|6b|) that prescribes £(t). It 
is not difficult to realize that to compute d£/dt at 2PN 
and 2.5PN orders, one only needs to know conservative 
orbital dynamics to 2PN order. However, computations 
that lead to d£/dt at 3PN and 3.5PN orders require con- 
servative orbital dynamics to 3PN order. We expect that 
the change in the conservative dynamics from 2PN to 
3PN order while computing d£/dt is the dominant rea- 
son for the observed gap in </>pn — </>nr between 2.5PN 
and 3PN orders in the bottom panel of Fig. [1] 

We provide Fig.[2]to argue that monotonic convergence 
to the 'exact' GW phase, exhibited by the TaylorEt tem- 
plates is rather insensitive to small numerical eccentric- 
ities present in NR evolutions. We obtained Fig. fj] by 
employing NR puncture evolutions having an initial ec- 
centricity e ?s 0.008. The oscillatory nature of </>pn — <^nr 
clearly demonstrates that we are indeed dealing with NR 
simulations that possess some tiny orbital eccentricity. It 
is also evident from Fig. [5] that 0pn — 0nr at any given 
PN order is not that substantially different from what 
is given in the bottom panel of Fig. [2j the accumulated 
phase difference was A</> = —1.18 radians for the data 
with negligible eccentricity, and is Acf> = —1.16 for the 
more eccentric data. 
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FIG. 1: Accumulated phase disagreement between NR and 
PN evolutions for TaylorTl, T4 and Et approximants at 2PN, 
2.5PN, 3PN and 3.5PN reactive orders. 
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FIG. 2: Accumulated phase disagreement between NR results 
and the PN approximant TaylorEt, as in the lower panel of 
Figure [l] with the difference that we now use numerical data 
from a simulation with slightly larger eccentricity, e ~ 0.008. 



Let us explore another unique feature of TaylorEt ap- 
proximant relevant for making contact with NR-based 
black-hole binary evolutions. A close inspection of 
Eqs. (O defining TaylorEt GW phase evolution reveals 
that if one can provide bounding values for £, it is possi- 
ble to obtain the associated <f>(t) at various PN orders. 

One can estimate the binding energy on one time slice 
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FIG. 3: The same comparison as in Figure [2j but we 
now match the phase and binding energy such that Eb — 
— 0.01383M at the matching time. 

of a numerical black-hole spacetime by calculating the 
difference between the total energy of the spacetime, the 
Arnowitt-Deser-Misner (ADM) mass, and the total mass 
of the two black holes. This estimate of the binding en- 
ergy was introduced in [29| . and is usually performed on 
the initial time slice, and as such includes the energy in 
the "junk" radiation associated with the standard initial- 
data choices. However, the junk radiation quickly radi- 
ates away, and then we may make a potentially more 
reliable estimate of the binding energy by calculating 

Eb = Eadm — Mi — M 2 — E ra d, (7) 

assuming that one has accurately measured the energy 
lost in gravitational-radiation emission, E ra d- We find 
that we can estimate the radiated energy with sufficient 
accuracy to place an uncertainty on the binding energy 
of less than 2%. 

When we perform a matching with respect to bind- 
ing energy, such that the binding energy and phase are 
equal at some time (in Figure [3J E\> = — 0.01383M for 
both NR and PN waveforms at the matching time), the 
GW frequencies of the NR and PN waveforms are not the 
same at that time. Therefore the time derivative of the 
phase disagreement is non-zero, and the phase disagree- 
ment immediately grows linearly as we move backwards 
in time. 

The curves in Figure [3] possess a turning point that 
makes an estimate of the accumulated phase error am- 
biguous. However, we should point out that this ambi- 
guity exists in all NR-PN phase comparisons, and this is 
clear in the figures shown in Refs. Q and ■ It is always 
possible to perform a matching at a different time (and 
therefore a different frequency or, in this case, binding 
energy), and to find a different accumulated phase dis- 
agreement, depending on the new location of the turning 
point. The comparison between the NR and PN phases 
is cleanest when a matching point is chosen such that 
there is no turning point, as in the analysis presented 
earlier and demonstrated in Figure [TJ Despite these am- 
biguities, we once again see only monotonic changes in 
the accumulated phase error when comparing NR results 
with the TaylorEt approximant. 



FIG. 4: Accumulated phase disagreement between NR and 
PN results, for each of the three approximants, TaylorTl, 
TaylorT4 and TaylorEt at 2PN, 2.5PN, 3PN, 3.5PN orders. 
At 2.5PN the TaylorTl and TaylorEt points are on top of each 
other. Note that the disagreement between NR and TaylorEt 
decreases monotonically as the reactive PN order is increased. 



V. DISCUSSION 

In this paper, we have compared GW phase evolu- 
tions associated with non-spinning equal-mass NR based 
black-hole binary inspiral lasting nine orbits with three 
different prescriptions based on the TaylorTl, T4 and 
Et approximants. We verified that GW phase evolution 
prescribed by TaylorT4 at 3.5PN order agrees very well 
with its NR counterpart. Further, </>pn — 0nr associated 
with TaylorTl and T4 approximants fluctuates as we in- 
crease PN order from 2PN to 3.5PN in steps of 0.5PN. 
However, the recently introduced TaylorEt approximant 
displays an intriguing feature, namely, 0pn — 0nr de- 
creases in a monotonic fashion as we increase the PN 
order responsible for reactive GW phase evolution. We 
also infer that 0pn — 4>nr associated with the TaylorEt 
approximant at 3.5PN order should be tolerable for GW 
data analysis purposes. Figure 2] summarizes our NR 
versus PN comparisons. 

The monotonic convergence to the fully general- 
relativistic GW phase evolution, exhibited by the Tay- 
lorEt templates, along with its other attributes listed 
earlier, makes it an interesting candidate for GW tem- 
plates that model inspiral, or for producing hybrid PN- 
NR waveforms as in [3(1 El, H3 ■ 

The observation that TaylorEt templates should be 
more efficient in capturing GWs from compact binaries 
having small orbital eccentricities in comparison with 
TaylorTl and T4, demonstrated in Ref. [lj] while re- 
stricting the radiation reaction to dominant contribu- 
tions, makes it a promising candidate to search for GWs 
using data from GEO-LIGO and VIRGO. Data analysis 
implications of GW polarizations h x ,+(t), evolving un- 
der the TaylorEt prescription, relevant for ground and 
space-based GW interferometers arc currently under in- 
vestigation. 

In future work we plan to make similar comparisons 
involving unequal-mass and spinning binaries. 
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